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Monolayer M0S2: trigonal warping, 'T- valley" and spin-orbit coupling effects 
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We use a combined ab-initio calculations and k ■ p theory based approach to derive a low-energy 
effective Hamiltonian for monolayer M0S2 at the K point of the Brillouin zone. It captures the 
features which are present in first-principles calculations but not explained by the theory of Xiao et 
al. [Phys Rev Lett 108, 196802 (2012)], namely the trigonal warping of the valence and conduction 
bands, the electron-hole symmetry breaking, and the spin-splitting of the conduction band. We also 
consider other points in the Brillouin zone which might be important for transport properties. Our 
findings lead to a more quantitative understanding of the properties of this material in the ballistic 
limit. 
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Introduction. Transition metal dichalcogenides are 
emerging as promising new materials for applications in 
electronics and optoelectronics 1 . In particular, mono- 
layer molybdenum disulphide (M0S2) has recently re- 
ceived significant attention experimentally 2 - - — as well as 
theoretically-^ 3 - - — . It may become the material of choice 
for field-effect transistors with high on-offratic 2 . In addi- 
tion, the strong spin-orbit coupling, the coupling between 
the spin and valley degrees of freedo m 13 ' 15 , and their 
effect on the exciton photoluminescense have sparked 
strong interest! 3 - - —. In light of the growing interest to 
this material, an accurate yet reasonably simple model 
that describes its band structure and electronic proper- 
ties is highly desirable. 

Following the important work in Refs. [12| and [ll| , 
Xiao et al*£ have recently introduced a tight-binding 
model which assumes that at the K point of the Brillouin 
zone (BZ) it is sufficient to take into account the d z 2 (for 
the conduction band) and d xy , d x i_ y i (for the valence 
band) atomic orbitals of the Mo atoms. Neglecting the 
spin-orbit coupling (SOC), Xiao et al. found an effective 
Hamiltonian of the form Hq = Hvo(Tk x (T x + k y a y ) + %<7 Z 
where cr x ,y, z denote Pauli matrices, A is the energy gap, 
vq plays the role of "Fermi- velocity" and r = 1 ( — 1) 
for valley K (K'). H$ describes massive Dirac parti- 
cles, in other words, it is a monolayer graphene Hamil- 
tonian with a staggered sublattice potential. While it 
seems to explain many experimental observations at least 
qualitatively- 5 - - — , certain limitations of this model can al- 
ready be appreciated by looking at Figs. QJb) and (c). 
The dispersion predicted by H is isotropic and possesses 
electron-hole symmetry regarding the valence and con- 
duction bands. As one can see in Figs. QJb) an d (c), 
which show the results of first-principles calculations, ex- 
cept in the immediate vicinity of the K point the disper- 
sion is not isotropic: a trigonal warping (TW) of the 
isoenergy contours can clearly be seen. In comparison to 
monolayer graphene we note that its low-energy disper- 
sion is isotropic on the energy scale of 1 eV, whereas in 



MoS 2 the TW is already observable at w 0.08 eV below 
the valence-band edge. Furthermore, ab initio calcula- 
tions predict different effective masses for electrons and 
hole o 18 i 19 , which obviously breaks the electron-hole sym- 
metry. We also note that both our density functional 
(DFT) calculations in the local spin density approxima- 
tion (LSDA) and the GW computations of Ref. [18] in- 
dicate that there is a relatively small (compared to the 
corresponding splitting in the valence band) but finite 
spin-splitting of about 3-4 meV in the conduction band 
at the K point, which cannot be explained in the theoret- 
ical framework of Ref. [13]. Finally, as one can observe 
in Fig. QJa) (see also the calculations of Refs. |18l - r20l ]) 
the valence-band maximum (VBM) at the V point is ac- 
tually very close in energy to the valence-band edge at 
the K point. While the exact value of the band maxi- 
mum seems to depend on the particular computational 
method that is used (c.f. Fig. 1(a) in Ref. [131 an d Fig. 3 
in Ref. [19J]), it is clear that at finite temperatures in 
hole-doped samples states at both K and T points will 
participate in transport (for other points of interest in 
the BZ see Appendix [Cjl . These observations call for a 
more exact model for the band structure of M0S2 . 

Using our ab initio computations and the kp theory 3 ^ 
we show that the observed TW of the bands can be de- 
scribed by a four-band generalized bilayer-graphene-type 
Hamiltonian and the TW is understood as a consequence 
of the coupling of the valence band (VB) and conduction 
band (CB) to other (remote) bands. From the four-band 
model we derive an effective two-band model and ob- 
tain the parameters that enter the model from fitting to 
our DFT computations. SOC also plays an important 
role in the low-energy physics of M0S2, but the SOC 
Hamiltonian cannot be obtained from results on bilayer 
graphen e 29 ' 30 . We find that a consistent description re- 
quires a seven-band model (or fourteen-band model, in- 
cluding the spin) at the K point and a six-band model at 
the r point. We derive an effective low-energy Hamilto- 
nian which takes into account the effects of the SOC, in- 




FIG. 1: a) Spin-resolved band structure of M0S2 from DFT 
LSDA calculations, b) Contour plot showing the isoenergy 
contours of the valence band (for zero SOC) from DFT cal- 
culations at the K point of the BZ (symbols) and as obtained 
from Eq. @ (solid lines), c) the same as in b) for the con- 
duction band, ao is the lattice constant. 



eluding the spin-splitting of the CB, which, to our knowl- 
edge has not been considered before. We also discuss 
how the effective mass and various other parameters of 
the model depend on the SOC. 

We start with the derivation of the spinless effective 
Hamiltonian, i.e. we neglect the SOC. This is already 
sufficient to explain the TW of the isoenergy contours. 
We make repeated use of various symmetries of the crys- 
tal structure, the two most important of which for our 
purposes are the rotational symmetry by 2tt/3 around 
an axis perpendicular to the plane of M0S2 (we denote it 
by C3) and the horizontal mirror plane ah- (For the full 
point group symmetry see Appendix fXJ) The derivation 
of the effective Hamiltonian relies on our DFT calcula- 
tions, which, in addition to the band structure, provide 
us with the projection of the DFT wave functions onto 
atomic orbitals at high symmetry points of the BZ. This 
helps us to identify the symmetries of the bands, which 
is necessary to obtain the general form of the effective 
Hamiltonian. As an example we consider the valence 
and the conduction bands at the K point of the BZ. 
Similarly to Refs. 01113 



we find that here the VB is 
predominantly composed of the d x 2_ y 2 and d xy atomic 
orbitals centered on the Mo atoms, which are symmetric 
with respect to ah- Since the VB is non-degenerate at 
the r point, compatibility relations require that at the 
K point it transforms as the A' representation of the 
group Csh, which is the small group of the wave vector 
at K. We denote the wavefunction of the VB by 1^^/), 
and hereafter we use the notation |^*) for the wavefunc- 
tions of various bands, where b denotes the band and 



\i the pertinent irreducible representation (irrep). The 
CB is predominantly composed of d z 2 orbitals of the Mo 
atom o 13 i 19 ~ — , which means that the CB wavefunction 
|\&£, ) is also symmetric with respect to 07, and trans- 
forms as the E[ irrep of C^h- Similar considerations al- 
low us to obtain the symmetries of all bands at the K 
point, even when different orbitals from different atoms 
are admixed (see Appendix |A"]) . In a minimal model that 
captures TW, in addition to the VB and the CB, there 
are two other important bands, both of which are even 
with respect to ah'- the second one above the CB, which 
we denote by CB+2 and whose wavefunction is \^ C j^ 2 ) 
and the third one below the valence band (VB-3) with 
wavefunction I'F^T 3 ). The other bands between the VB- 
3 and CB+2 are antisymmetric with respect to the mir- 
ror plane of M0S2 and therefore they do not couple to 
the VB and the CB. In k • p theory the Hamiltonian 
Hk-p = — k • p is considered a perturbation (m e is the 
bare electron mass) and one uses first-order perturbation 



theory in the basis of {|v&^/), |* 
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Using the notation k = q + K, the perturbation can 
be rewritten as % k . p = \-^{q+P- + Q-P+) = "H kp + 
% k . p , where the operators p± are defined as p± = 
p x ± ipy and similarly q± = q x ± iq y . The matrix el- 
ements of "Hkp are constrained by the symmetries of 
the system. For example, considering the rotation C3, 
the relation (^ v A ,\p + \ty 
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= <^,|C3 t C3P + C 3 t C 3 |^t 2 > 



should hold. Since (^\'\ c z 
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that (n'l^.pint 2 ) = ^ i47r/3 (n>\K. P mf)> which 

means that this matrix element must vanish. By con- 
trast, (^A'l^kpl'^'ir 2 ) = 74 1S finite and one can prove 
that it is a real number. Similar considerations can be 
used to calculate all matrix elements. Finally, in the basis 
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at the K point is given by 
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^ C J, )} the k • p Hamiltonian 



#kp = 



/ £v 739- 729+ 749+ 

739+ £c 759- 769- 

729- 759+ e»>-3 

V 749- 769+ e c+2 



(1) 



where 7$ are the matrix elements of Hk p in the above 
mentioned basis and e v ,£ c ,e v —3,e c +2 are band-edge en- 
ergies. The matrix elements between I^^T 3 ) and |\fr^> ) 
vanish due to symmetry. We note in passing that the 
Hamiltonian ([T]) can be considered a generalized bilayer 
graphene (BLG) Hamiltonian. This can be seen by ro- 
tating the well known bilayer graphene Hamiltonian into 
a basis where the basis functions transform according to 
the irreps of the small group of BLG (see Appendix lB|l . 
The k • p Hamiltonian at the K' point can simply be 
obtained by re-defining q± as q± = q x ± ir q y , where 
t = 1 (-1) for the K (K') point. As in the case of BLG, 
it is convenient to derive a low-energy effective Hamilto- 
nian from |T]) , which does not contain the contribution of 



the bands far from the Fermi energy. Using the Lowdin 
partitioning 35 we find 



Ho + H a 



H c s = H + H 3w + H as + H cu b, (2a) 
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(q-) 2 
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where the parameters a, /3, 7 and r\ can be expressed 
in terms of 7$ as defined in Eq. ((T|) and the band edge 
energies. The parameters a and /3 describe the breaking 
of the electron- hole symmetry that is apparent comparing 
Figs. [T](b) and (c), whereas 7 is responsible for the TW of 
the energy contours. The cubic term, H cuo is important 
to achieve a quantitative fit to the VB away from the 
K point. We note that symmetries allow for another, 
diagonal Hamiltonian which is ~ g 3 , but its effect has 
been found to be negligible. 

We used the VASP code 31 and the local density ap- 
proximation (LDA) of DFT to calculate the band struc- 
ture and isoenergy contours shown in Fig.[TJ To compare 
our k • p theory with the DFT calculations one has to de- 
termine the matrix elements 7$. These matrix elements, 
in principle, can also be calculated from DFT 27 . We 
found however, that the isoenergy lines calculated with 
parameters jf FT obtained from the numerical evaluation 
of (^ulPil^^/) using Kohn-Sham orbitals 3 - 2 - give a satis- 
factory agreement with DFT band structure only in the 
close vicinity of the K point. Therefore we used these 
7j DFT values as a starting point for a fitting procedure 
whereby we fitted the eigenvalues of the Hamiltonian @ 
[given in terms of 73, a, /3, 7, r\ and the band edge ener- 
gies e v , and e c ] along the TKM line to the CB and VB 
obtained from the DFT calculations. From the fitting we 
found the parameters 73 = 3.82 eVA, a = 1.72 eVA 2 , j3 — 
-0.13eVA 2 , 7 = -1.02 eVA 2 and 77 = 8.53 eV A 3 . The 
isoenergy contours calculated using these parameters [the 
solid lines in Figs. HJb) and (c)] capture well the TW of 
the band dispersion (c.f the results of the DFT computa- 
tions given by symbols) , which is more pronounced in the 
VB than in the CB. The agreement between the DFT re- 
sults and the predictions based on ^ is very good up to 
energies 0.16 eV below the VB maximum and above the 
CB minimum; for other energies the agreement is qualita- 
tive. The effective masses (along the TKM line) that can 
be inferred from these parameters are m^ « — 0.62m e 
for the VB and m 



An important feature of the band structure of M0S2 
which has received little attention so far is that the top 
of the valence band at the T point is very close in energy 
to the VBM at the K point±2r.2&; see also Fig (Ha). This 
means that for the VB transport properties the states 
close to the T point can also be important. The Hamil- 
tonian of this 'T-valley" can also be derived using the 
kp theory, along similar lines to the case of the K point. 
Note, however, that the group of the wave vector at the 
r point is D^h- Our DFT calculations show that here 
the VB is mainly composed of the d z i and s orbitals of 
the Mo atoms and the p z orbitals of the S atoms, which 
means that it belongs to the A x irrep of D3/,. The VB is 
coupled to the VB-3 and CB+1 bands which are doubly 
degenerate at the Y point. There is no coupling between 
the VB and the CB at the Y point: due to band cross- 
ings along the T — K line the CB becomes antisymmetric 
with respect to ah- Upon performing the Lowdin parti- 
tioning we find that the dispersion is isotropic and can be 
well described by Hr = " r , where the effective mass 

m^g = — 3.647m e is found by fitting the band structure, 
which is in good agreement with Ref. [22j. The impor- 
tance and role of the Y point in the transport properties 
of the VB would require further experimental work. We 
expect that in disordered samples due to their large ef- 
fective mass and hence low mobility the contribution of 
these states to the transport is small, but they may be 
important in ballistic samples and in quantum-Hall mea- 
surements. 

The description of the system becomes more compli- 
cated if one takes into account the SOC as well. In the 
atomic approximation the SOC is given by the Hamilto- 
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0.48m e for the CB, which are 
in good agreement with the results of other DFT LDA 
calculation o 17 i 22 . Interestingly, we have checked by nu- 
merical calculations that although the TW is quite strong 
in the VB, its effect on the Landau levels is actually very 
small, and they can be calculated by neglecting H 3w in 

Eq. © .... " , „ „ ..,. , 

Berry-curvature and hence various Hall-conductivities^. show that (s, *^,|L_5 , + |*^„ , s) = A( VV ^ 1 - ) S + where 



Nevertheless, H 3w and H cub should affect the the fact that c 3 L ± C 3 



Here V(r) is the spherically symmetric atomic po- 
tential, L is the angular momentum operator and 
S = (S x ,S y ,S z ) is a vector of spin Pauli ma- 
trices S x , S y (with eigenvalues ±1). Note that 
L- S = L Z S Z + L + S- + L-S+, where L± = L x ± iL y 
and S± = 2^ x -*- iSy)' Let us introduce the spinful 
symmetry basis functions by |\^, s) — \^ b n) ® \s), 
where s = {t,i} denotes the spin degree of freedom, 
and consider first the K point of the BZ. Since L± 
transforms as the E" irrep of C^, there can be non- 
vanishing matrix elements of H^ between states that are 
even/odd with respect to ah- Therefore we considered a 
seven-band model (without spin) which contains every 
band between VB-3 and CB+2, i.e., we consider the basis 

l^ ^ 2 , s}} . The symmetries ah and C3 of the system 

here also help us to find the non-zero matrix ele- 
ments of H^. For example, one can make use of 



e Tl2 ^/ z L± and therefore 



A(v,v—i) i s a constant, whereas (s, ^\,\L + S-\'9 V „ 1 , s) — where our DFT computations give a spin-splitting of 



(s, ^ V A ,\L Z S Z \^ V ,, , s) = 0. By calculating the matrix 

2 
iJ°o in the above mentioned basis and diagonalizing 

the Hamiltonian Hd + H^ where Hd is a diagonal 
matrix containing the band-edge energies, one obtains 
the eigenstates \^ b ,, s), which in general turn out 
to be linear combinations of a symmetric | "J^ , s) and 
an antisymmetric |\IJ*, , s) wavefunction with different 
weights. In our notation the new eigenstates l^.^s) 
inherit the band index b and spin index s from the 
state whose weight is larger in the linear combination 
that makes up \ty b „/,s). This assignment of the band 
index and spin quantum number is possible because 
the typical energy scale of the SOC (the upper limit of 
which is the splitting of the valence band sw 145meV, 
sec below) is smaller than the typical band separation, 
i.e., the bands are not strongly hybridized by the SOC. 
The diagonalization of the Hamiltonian can be done 
analytically in the approximation that couplings of up 
to next-nearest-neighbour bands are kept and more 
remote couplings, e.g., between |^"7 3 ,s) and \^f c X ,s) 

E 2 A 

are neglected. All eigenstates are non-degenerate, 
as expected, since the double group of C^h has only 
one-dimensional representations. With the new eigen- 
states | *f? b , , s) one can repeat the k • p calculation, 

and since |$' /,s) is an admixture of symmetric and 
antisymmetric states, there will be more non-zero matrix 
elements of the Hk P Hamiltonian than there were in 
the case of zero spin-orbit coupling; see Eq.([T|). Finally, 
using the Lowdin partitioning we can derive an effective 
low-energy Hamiltonian for the spinful valence and 
conduction bands. Since the calculations are quite 
lengthy, we only give the most important results here 
and concentrate on the zero-magnetic field case. The 
Landau-level problem in the presence of SOC and the 
related question of the effective g-factor of monolayer 
M0S2 will be discussed elsewhere. 

We will work in the basis of {|#" ^/.t), 1*1, „»4 
),|*°; c ,«,t)j|*i' a" A)} and start with the diagonal 

E 2 ,E 1 E 2 ,A 

and q independent part of the SOC Hamiltonian, i.e., 
we consider the SOC dependence of the band edge ener- 
gies. According to our kp calculations, the spin splitting 
of the VB is given by H% v « -tA\S z + ^(1 + tS z ), 
where our DFT computations give a spin-orbit gap of 
2A? - A| « 146meV. The first term, H% v « -rAJS* 
was first obtained in Ref. [13l |. whereas the second term, 
which is expected to be much smaller, comes from the 
coupling of spin- up (spin-down) band to VB-1 at the K 
(K 1 ) point. (The coupling of the spin-down (spin-up) 
band of the VB to other bands is much weaker.) The 
spin splitting of the CB, although noted in Ref. [lg|, has 
not yet been discussed in the literature. It can be de- 
scribed by the Hamiltonian 
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3meV. Although this is a 

small effect compared to the spin-splitting in the VB, 
spin-splittings of similar magnitude have recently been 
measured in, e.g., carbon nanotube quantum dots 2 ^. The 
spin-splitting in the CB originates from the SOC of the 
CB to the VB-1 and CB+1 bands, and it is a consequence 
of the hitherto neglected off-diagonal SOC terms, related 
to the ~ L-S++L+S- part of %**. Our results therefore 
show that the spin- valley coupling is present not only in 
the VB 13 but also in the CB. 



Regarding the effect of SOC on the q-dependent terms 
in Eq. ([2]), we find that bands with different spin in- 
dices, e.g., {|*' />B «,t) and 1*^,^,1) or 1*^,1}} 
do not couple to each other. After folding down the full 



seven-band k • p Hamiltonian, in the basis of {l^^, p « , t 

1*^' vol)} the effective 

a 2 ,A 



>,|^, x ,t>}and{|^, x 4>, 
Hamiltonian is still of the form of Eq. ([2]), but in general 
with different parameters 7J , o^- 1 , /3^^- 1 , ^^ and 
^t(-l-) f or thg spi n - U p (spin-down) bands. By fitting our 
SOC-resolved DFT calculations we find that 73 and 0, 
hence m^ are basically not affected by the SOC. The 
effective masses in the VB are slightly renormalized by 
the SOC, leading to m^g ~ 0.65m e and m v ^ « 0.5877i e , 
i.e., a difference of roughly 5% — 6% with respect to the 
zero SOC case. 



(4) 



Finally, we briefly discuss the effect of SOC on the 
states at the T point of the BZ. In contrast to the K 
point, here the valence band remains degenerate even if 
we take into account SOC (see Fig[T]). This can be un- 
derstood from general group theoretical arguments: the 
pertaining double group has two-dimensional, hence de- 
generate representations. The dispersion for each spin 
can be described by a parabolic dependence on k and we 
find that the effective mass is basically unchanged with 
respect to the zero SOC case. 

Conclusions. We have derived a low-energy effective 
Hamiltonian for monolayer M0S2 at the K point of the 
BZ, which takes into account effects that are present in 
first-principles calculations but have not hitherto been 
discussed. The TW and spin-splitting of the bands 
should be readily observable by spin and angle resolved 
photoemission spectroscopy. We have also considered the 
states at the T point of the BZ, which can be important 
for transport properties of hole-doped samples as well as 
for various scattering and relaxation processe s 23 ' 26 , be- 
cause scattering from the K to the T point does not re- 
quire a simultaneous valley- and spin-flip. 
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Note added. During the preparation of this manuscript 
two related preprints have appeare d 36 ' 37 , where some of 
the results that we present here have also been obtained. 



Appendix A: Character tables and basis functions 

In Fig. [5] we show a top view of the monolayer M0S2 
lattice. The pertinent point groups to understand the 
band structure of monolayer M0S2 are D 3 h and C^h . The 
former is the group of the wave vector at the T point, 
the latter at the K point of the Brillouin zone (BZ). 
The symmetry operations that generate these groups are 
three-fold rotation C3 around an axis perpendicular to 
the plane of M0S2, a horizontal mirror plane ah perpen- 
dicular to the three- fold axis and in the case of D^h, three 
two-fold rotation axis C 2 that lie in the horizontal mirror 
plane. 



^^.. 




symmetry properties, the analysis of the bands in terms 
of atomic orbitals, together with band compatibility re- 
lations, gives us the irreps that can be assigned to each 
band. 

As an example we consider the valence band (VB) at 
the K point. Here the VB is predominantly composed of 
the d x 2_ y 2 and d xy atomic orbitals centered on the Mo 
atoms, which are symmetric with respect to ah- The two 
Bloch functions that can be formed from these orbitals 
and which reflects the three-fold rotational symmetry are 
l*2%(k)) = ^E n e MR " +t ^Y 2 ±2 (r-(R n +t Mo )) 

where V 2 ±2 (r) — (d x 2_ y 2(v) ±id xy (r)) /\/2. The small 
group of the wave vector at the K point is C$h ■ One 
can show that ^^(K)) transforms as the E' 2 irrep of 
this group whereas ^^^(K)) transforms as the A' irrep. 
Since the VB is non-degenerate at the T point, compati- 
bility relations require that at the K point it transforms 
as the A' irrep. 

Table Q] shows the characters and irreps for D^h, perti- 
nent to the r point, while Table UT1 shows symmetry prop- 
erties of the rotating orbitals based Bloch wave functions 
and that which band they contribute to. The conduction 
band is denoted by CB, the valence by VB, the first band 
above the CB by CB+1, the first band below the VB by 
VB-1, and so on. 



FIG. 2: Top view of the M0S2 lattice. Mo atoms are indicated 
by grey (solid line) circles, S atoms by yellow (dotted line) cir- 
cles. The lattice vectors ai = ^(l,\/3) and a 2 = ^-(1, — s/3) 
are also shown (ao = 3.129A is the lattice constant). 

By projecting the plane-wave basis used in our DFT 
computations onto atomic orbitals one can obtain the 
decomposition of each band in terms of atomic orbitals 
Q r l, where 77 = {Mo, S1,S2} denotes whether the given 
orbital is centered on molybdenum (Mo) or on one of the 
sulphur (SI, S2) atoms in the unit cell, and the lower in- 
dex v = {s,p Xl p yi p z ,d z 2 1 d xy ,d XZl d xy ,d x 2_ y 2} indicates 
the type of orbital. To take into account the three-fold 
rotational symmetry of the system, one should use linear 
combinations of these orbitals to form the rotating or- 
bitals Y™ , which are proportional to spherical harmon- 
ics. We then consider the transformation properties of 
the Bloch wave functions formed with the rotating or- 
bitals: 



l^, m (k,r)) = 



1 



J2e lk < R " +t ^ Yr (r-{R n +t n }). 



Here the summation runs over all lattice vectors R n and 
t Mo = f (1, ~^) T > tsi = t S2 - f (1, ^) T give the po- 
sition of the Mo and S atoms in the (two-dimensional) 
unit cell with ao = 3. 129 A being the lattice constant (see 
also Figf2]) and k is measured from the T point of the 
BZ. We then identify the irreducible representations (ir- 
reps) according to which ^^(k, r)) transform at the 
high symmetry points T and K of the BZ. Since hy- 
bridization between different orbitals will preserve the 



6m2 (D 3h ) 


E on 2C3 25*3 3C2 3<r„ 


A[ 


111111 


A' 2 


11 1 1-1-1 


A'{ 


1-11-11 1 


A'i 


1-11 -1 -1 -1 


E' 


2 2-1-10 


E" 


2-2-11 



TABLE I: Character table and irreps of the group 6m2 (D3 



irrep 



X 



A'i 



E 1 



E" 



basis functions 



*o M o°), |*£o°>, -75 Cl«-f,o> - l*f,o» 



Ko>, £ (!*&> - l*i.o» 



{|*2f2 >,|*2^2>} 

{71 (I*m> - l*f,i» - ^ (l*f,-i) - l*f 2 -i»} 



{|*ft">,|*a£i» 

{72 0*m> - 1*0)) , 75 (l*i'.-i) - l*f 2 -i»} 



band 



VB 



VB-2 



VB-3 



VB-1 



TABLE II: Basis functions for the irreps of the small group 
Dzh of the r point. {. . . } denote the partners of the two- 
dimensional representations. The rightmost column shows 
that to which band the basis functions contribute. 



Table IIIII gives the characters and irreps of C^h > perti- 
nent to the K point and Table IIVI the rotating orbitals 
based Bloch wave function which transform as the irreps 
oiC 3h . 



G(C sh ) 


E C 3 Ct 


Oh 


Sz 


°~hC 3 


A' 


1 1 1 


1 


1 


1 


A" 


1 1 1 


-1 


-1 


-1 


E{ 


1 2 


1 


U) 
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UJ 


E' 2 


1 W UJ 


1 


u? 


UJ 


E'{ 


1 UJ UJ 


-1 


— UJ 


-UJ 2 


E'i 


1 2 

1 UJ UJ 


-1 


2 

— UJ 


— UJ 



TABLE III: Character table of the group 6 (C 3 



irrep 


basis functions 


band 


A' 


|*tf°2>, ^ (l^f, 1 -!) + l^f, 2 -!)) 


VB 


A" 


I*ft°>, 73 (l*!, 1 -!)-!*!, 2 -!)) 


CB+1 


E[ 


1*^°),^ (|*f,i> + l*f, 2 i» 


CB 


E' 2 


l*&°>, ^(I*f.o)-I*f 2 o)) 


VB-3 
CB+2 


E'l 


lO^(KD-l^)) 


VB-2 


E'i 


l*r-i),7f(l*f,o> + l*f,o)) 


VB-1 



TABLE IV: Basis functions for the irreducible representations 
of the small group of the K point. The rightmost column 
shows that to which band the basis functions contribute. The 
basis functions for the K' point can be obtained by complex- 
conjugation. 



hopping between the atoms Bl and A2. One can per- 
form a unitary transformation which rotates the Hamil- 
tonian JBI} into the basis {\$ai), \®A2), |*Ei), \^e 2 )} 
where the basis functions 1^^) transform as the irreps 
fi = {Al, A2, E} of the small group of the K point, which 
is 7?3 in this case. One finds 



/ 



H, 



BG 
kp 



V 



~7i 





-V04P+ 


voiP- \ 





7i 


V04P+ 


VQ4P- 


-V04.P- 


V 4P- 





V3P+ 


VQ4P+ 


VQ4P+ 


V3P- 


o / 



(B2) 



where v oi = 4=(v -V4,), and ^04 = -^{v + v 4 ). This 
Hamiltonian is characterized by the three hoppings V3, 
W04, and 1T04, and three band edge energies —71, 71, 
and (degenerate). The well known low-energy effec- 
tive Hamiltonian of bilayer graphene^ can be obtained 
by projecting out the states \^ai) and | v E'a2)- 

The Hamiltonian of monolayer M0S2 has the same 
structure as (|B2[) but is characterized by five different 
hoppings and four different band-edge energies; in this 
sense it is a generalization of (|B2|) . 



Appendix C: Q-point minimum in the conduction 
band 



In this section we briefly discuss whether it is impor- 
tant to consider the minimum at the Q point in the con- 
duction band (see Fig. [3]). 




Appendix B: Bilayer graphene Hamiltonian 

The k • p Hamiltonian of bilayer graphene^ at the K 
point of the BZ, in the basis of {A2, Bl, Al, B2} sites is 
given by 



ttBG _ 

H k P - 



/ 71 v 4 p + v p- \ 

71 V p + V4P- 

vaP- «oP- v 3 p + 

\ v p + v 4 p + v 3 p^ / 



(Bl) 



where we have chosen the on-site energies to be zero, 
P± = Px ± ipy, the velocities vq , V3, and V4 depend 
on intra- and interlayer hoppings, and 71 is the direct 



FIG. 3: Conduction band of M0S2 from DFT calculations 
using the HSE06 functional (red, solid line) and the LDA 
(green, dashed). 

To this end we compare the conduction band calcu- 
lated in local density approximation (LDA) and using 
the HSE06 hybrid functional 39 . The latter, while being 
computationally much more demanding than the LDA, 
has been shown to improve the accuracy of electronic 
structure calculations for many semiconductors^. 

As one can see in Fig. [3J there are two main differ- 
ences between the results of calculated with the HSE06 
functional (red, solid line) and LDA (green, dashed 
line). Firstly, there is an up-shift of the HSE06 con- 
duction band with respect to the LDA one, leading to 



a larger band gap at K. Secondly, the minimum at the 
Q point is much higher in energy (and becomes shal- 
lower) than the minimum at K in the case of HSE06 
calculations. In particular, the difference between the 

0.405 eV for HSE06 and 



minima is E. 



E, 



hybrid 



E 



hybrid 
K 



LDA 



E\f A = 0.09 eV for LDA. (Note, that the 
LDA calculations of Ref. [l7[ give ps 0.2 eV). For com- 
parison, the difference between the valence band max- 
ima are £^ brid - £ r iybrid = 0.058 eV for HSE06 and 
E)p A - Ep> A = 0.12 eV for LDA. Therefore, regard- 
ing transport properties, for p-doped samples states at 



the F point are more important than the states at Q for 
the n-doped case. 

We note that both the increase of the band gap at the 
K point and the up-shift of the minimum at the Q point 
are in qualitative agreement with the GW calculations of 
Ref. [lj| • The importance of the Q point minimum can 
hopefully be determined when more accurate measure- 
ments of mobility become available, because the phonon- 
limited mobility depends quite sensitively on the energy 
separation of Ek — Eq (for details see Ref. |23|). 



* e-mail: andor.kormanyos@uni-konstanz.de 

1 Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, 
and M. S. Strano, Nature Nanotechnology 7, 699 (2012). 

2 B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, 
and A. Kis, Nature Nanotechnology 6, 147 (2011). 

3 A. Splendiani, L. Sun, Y. Zhang, T. Li, J. Kim, Ch.-Y. 
Chim, G. Galli, and Feng Wang, Nano Letters 10, 1271 
(2010). 

4 G. Eda, H. Yamaguchi, D. Voiry, T. Fujita, M. Chen, and 
M. Chhowalla, Nano Letters 11, 5111 (2011). 

5 H. Zeng, J. Dai, W. Yao, D. Xiao, and X. Cui, Nature 
Nanotechnology 7, 490 (2012). 

6 K. F. Mak, K. He, J. Shan and T. F. Heinz, Nature Nan- 
otechnology 7, 494 (2012). 

7 T. Cao, G. Wang, W. Han, H. Ye, Ch. Zhu, J. Shi, Q. Niu, 
P. Tan, E. Wang, B. Liu and J. Feng, Nature Communi- 
cations 3, 887 (2012). 

8 S. Wu, J. S. Ross, G. Aivazian, A. Jones, Z. Fei, G.-B. 
Liu, W. Zhu, P . Xiao, W. Yao, D. Cobden, and X. Xu, 
arXiv: 1208.6069 (unpublished). 

9 G. Sallen, L. Bouet, X. Marie, G. Wang, C. R. Zhu, W. 
P. Han, Y. Lu, P. H. Tan, T. Amand, B. L. Liu, and B. 
Urbaszek, Phys. Rev. B 86, 081301 (2012). 

10 M. Buscema, M. Barkelid, V. Zwiller, H. S. J. van der Zant, 
G. A. Steele, and A. Castellanos-Gomez, Nano Lett. 13, 
358 (2013). 

11 L. F. Mattheis, Phys. Rev. B 8, 3719 (1973); 

12 R. A. Bromley, R. B. Murray, and A. D. Yoffe, J. Phys. C: 
Solid State Phys. 5, 759 (1972). 

13 D. Xiao, G.-B. Liu, W. Feng, X. Xu, and W. Yao, Phys. 
Rev. Lett. 108, 196802 (2012). 

14 H.-Zh. Lu, W. Yao, D. Xiao, and Sh.-Q. Shen, Phys. Rev. 
Lett. 110, 016806 (2013). 

15 X. Li, F. Zhang, and Q. Niu, Phys. Rev. Lett. 110, 066803 
(2013). 

16 Z. Y. Zhu, Y. C. Cheng, and U. Schwingenschlogl, Phys. 
Rev. B 84, 153402 (2011). 

17 K. Kaasbjerg, K. S. Thygesen, and K. W. Jacobsen Phys. 
Rev. B 85, 115317 (2012). 

18 T. Cheiwchanchamnangij and W. R. L. Lambrecht, Phys. 
Rev. B 85, 205302 (2012). 

19 H. Shi, H. Pan , Y.-W. Zhang, and B. I. Yakobson 
larXiv:1211.5653l (unpublished) . 

20 J. K. Ellis, M. J. Lucero, and G. E. Scuseria, Appl. Phys. 
Lett. 99, 261908 (2011). 

21 E. S. Kadantsev and P. Hawrylak, Solid State Communi- 



cations 152 909 (2012). 

22 W. S. Yun, S. W. Han, S. Ch. Hong, I. G. Kim, and J. D. 
Lee, Phys. Rev. B 85, 033305 (2012). 

23 X. Li, J. T. Mullen, Zh. Jin, K. M. Borysenko, M. B. 
Nardelli, and K. W. Kim. larXiv:1301. 77091 (unpublished') . 

24 K. Dolui, I. Rungger, and Stefano Sanvito, 
larXiv:1301. 24911 (unpublished) . 

25 R. Roldan , E. Cappelluti, and F. Guinea. larXiv: 1301 .48611 
(unpublished). 

26 Yang Song and Hanan Dery, arXiv:1302. 36271 (unpub- 
lished) . 

27 Ch.-X. Liu, X.-L. Qi, H. Zhang, X. Dai, Zh. Fang and 
Sh.-Ch. Zhang, Phys. Rev. B 82, 045122 (2010). 

28 T. S. Jespersen, K. G. Rasmussen, J. Paaske, K. Muraki, 
T. Fujisawa, J. Nygard, and K. Flensberg, Nat. Phys. 7, 
348 (2011). 

29 A. Kormanyos and G. Burkard, Phys. Rev. B 87, 045419 
(2013). 

30 S. Konschuch, M. Gmitra, D. Kochan, and J. Fabian, Phys. 
Rev. B 85, 115423 (2012). 

31 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

32 The matrix elements of the momentum operator were eval- 
uated within the LDA using the CASTEP code^ 3 .. We used 
ultrasoft pseudopotentials, a plane- wave cutoff energy of 80 
Ha an artificial periodicity of 30 Bohr in the vertical di- 
rection and a 21x21 Monkhorst-Pack mesh. The optimised 
lattice parameter was 5.89 Bohr and the distance between 
the two layers of S atoms was 5.86 Bohr. 

33 S.J. Clark et al, Z. Kristallogr. 220, 567 (2005). 

34 M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group 
Theory, Springer- Verlag Berlin Heidelberg (2008). 

35 R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems, Springer- Verlag 
Berlin Heidelberg (2003). 

36 H. Rostami, A. G. Moghaddam, and R. Asgari, 
larXiv: 1302.59011 (unpublished). 

37 H. Ochoa, and R. Roldan. [arXiv:1303 . 5806 (unpublished). 

38 E. McCann and V. I. Fal'ko, Phys Rev Lett 96, 086805 
(2006). 

39 J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 
118, 8207 (2003); ibid 114, 219906 (2006). 

40 B. G. Janesko, T. M. Henderson, and G. E. Scuseria, Phys. 
Chem. Chem. Phys. 11, 443 (2009). 



